Extracting joint weak values from two-dimensional spatial displacements 
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The joint weak value is a counterfactual quantity related to quantum correlations and quantum 
dynamics, which can be retrieved via weak measurements, as initiated by Aharonov and colleagues. 
In this Rapid Communication, we provide a full analytical extension of the method described by 
Puentes et al. [Phys. Rev. Lett. 109, 040401 (2012)], to extract the joint weak values of single- 
particle operators from two-dimensional spatial displacements of Laguerre-Gauss probe states, for 
the case of azimuthal index |/| > 1. This method has a statistical advantage over previous ones 
since information about the conjugate observable, i.e., the momentum displacement of the probe, is 
not required. Moreover, we demonstrate that under certain conditions, the joint weak value can be 
extracted directly from spatial displacements without any additional data processing. 
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As the well-known examples of the Einstein-Podolsky- 
Rosen paradox and the which-path information in the 
Young double-slit experiment illustrate, quantum dy- 
namics and quantum correlations have eluded our intu- 
itive understanding for a long time. To capture these fea- 
tures, it is important to access the expectation value of 
products of two observables since such joint observables 
can contain information about quantum correlations and 
quantum dynamics. However, in the standard technique, 
the strong measurement of joint observables is difficult 
to implement because of the requirement of a nonlinear 
Hamiltonian [l[ and the state reduction resulting from 
the measurement backaction. In this regard, the Joint 
weak value is a good candidate, as seen in Refs. [2nali 
since this quantity can be experimentally obtained with 
a tiny measurement backaction from the second-order ef- 
fects of the time evolution via a local, linear, and single- 
particle interaction Hamiltonian. 

Historically, the weak value was initiated by Aharonov, 
Albert, and Vaidman (AAV) @, inspired from the two- 
time formulation of the quantum mechanical system . 
This formulation is characterized by the pre- and post- 
selected states of the system. When a system initially 
prepared in state |V'i) is post-selected in state the 
weak value of the observable A is defined as 



(V'fl^lV'i) 



(1) 



which can lie outside the spectrum of A and can even 
take an imaginary number. Recently, there have been 
various experimental realizations using the weak value 
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within the foundations of quantum mechanics, for exam- 
ple, in the direct measurement of a wavefunction Q , the 
which-path measurement in the Young double slit exper- 
iment [9|, the confirmation of the Heisenberg-Ozawa un- 
certainty relationship [loj , the confirmation of the Hardy 
paradox [llj, the violation of the Leggett-Garg inequal- 
ity [l3| , and the measurement of a geometric phase in an 
interferometer [l3| . 

There are also various proposals for extracting the 
weak value from experimental data (see for reviews El- 
Il6j). Originally, AAV restricted their attention to the 
standard von Neumann paradigm with a weak interac- 
tion Hamiltonian of the form H = gA (E) Px , where 5 is a 
small coupling constant and Px is the momentum observ- 
able of the probe state conjugate to the position observ- 
able X with [X, Px] = ih. Moreover, they assumed that 
the probe state was initially prepared in a fundamental 
Gaussian mode. In this case, we can determine the real 
and imaginary parts of the weak value from the spatial 
and momentum displacements of the probe state [a, ■ 

Recently, Resch and Steinberg Q proposed a measure- 
ment technique for the weak value of the joint observ- 
able as alluded before. They employed a two-dimensional 
Gaussian probe state and a weak coupling of the sys- 
tem observable A (B) with x (y) dimension of the probe 
state. By performing a second order expansion in the 
two-dimensional displacement of the probe state, they 
showed that it is possible to extract the real part of the 
joint weak value (A_B)w from the second-order spatial 
displacement, under the assumption of commuting ob- 
servables [A, B] = 0. This procedure, however, needs the 
real and imaginary parts of single weak values (A)w and 
(-B)w to calculate the joint weak value. Thus, we have 
to obtain full information of the probe wavefunction, i.e., 
not only the spatial displacement but also the momen- 
tum displacement, by taking the Fourier transform of the 
probe wavefunction. 
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In this Rapid Communication, we provide a simple 
method for extracting the joint weak value only from the 
spatial displacement of the two-dimensional probe state. 
This method has a statistical advantage over previous 
methods since full information of the probe wavefunc- 
tion is not required. The key idea is to employ Laguerre- 
Gauss (LG) modes for probe states, as was initiated by 
Puentes et al. Q. We extend this idea and provide the 
full description of the weak measurement for higher order 
LG modes with radial index p = to extract the joint 
weak values. 

The LG modes are given as the natural solutions of 
the paraxial wave-equation and characterized by a 
radial index p and an azimuthal index I. The modes have 
zero intensity at its center and an annular intensity dis- 
tribution. The wavefront of the LG modes is composed 
oi \l\ intertwined helical wavcfronts, with a handedness 
given by the sign of It has been shown that each pho- 
ton in the LG modes carries a quantized intrinsic orbital 
angular momentum Ih, in addition to the spin-like an- 
gular momentum ±/i associated with circularly polarized 
waves [l^. The LG modes have been created using var- 
ious experimental setups, e.g., using spatial light modu- 
lators 1201 and using the reflection on a conically-shaped 
mirror [21[. Furthermore, the LG modes have many ap- 
plications, for example, for achieving high efficiency of 
optical tweezers p^ . for reducing the thermal noise in- 
side gravitational- wave interferometers [2311 a nd for gen- 
erating entanglement with high efhciency [2J]. The am- 
plitude distribution of the LG modes with radial index 
p = is given as 



=iV{.T + i-sgn(/)y}l'l 



exp 



4(t2 



, (2) 



where a is the variance in the case of / = 0, sgn(-) is 
the sign function, and N is the normalization constant. 
When I = Q, Eq. ^ corresponds to a fundamental Gaus- 
sian mode. In this case, the amplitude distribution is 
factorable in two directions, x and y. When \l\ > 0, how- 
ever, it is no longer factorable, and this is a key factor 
for retrieving the joint weak values. 

Consider a weak interaction between the LG probe 
state = J dxdy(j)i{x,y)\x,y) and an initial state \ipi) 
of the system for the joint weak measurement of the 
obscrvables A and B. The total input state is {"^i) = 
I'tpi) ® |0i}. For the measurement process, we use the 
standard von Neumann paradigm with the local, linear, 
and single-particle interaction Hamiltonian as follows: 



(3) 



H = gSit -to)[A®P^ +B®Py 



where a coupling constant g is sufficiently small, and Px 
and Py are the momentum observables of the probe con- 
jugate to two commuting position observables X and Y ^ 
respectively. Here we have taken the interaction to be 
impulsive at time t ^ to and the same coupling constant 
g between x and y directions for simplicity. 



After the interaction between the system and probe 
states, we post-select the system in state lipi), resulting 
in the probe state: 



|0f) = (-'/'fie 



-^ig{A<S)P^+B<g)Py)/h 



(4) 



We denote the expectation value of an observable M in 
the initial probe state |0i) and the final probe state 
as 

We can calculate the displacement (M)fi by expanding 
the time evolution operator in Eq. (|4]) up to the second 
order in the coupling constant g: 



(Af)fi = (M)f - (M)i 

^ ^ lm{MHi)i + |j (^{HlAIHi)i - Re{MH2)i 



(6) 



where 



(7) 



Hl=Px{A)^ + Py{B)^, 

H2 = 4'(i2)w + Py{B^)^ + PxPy{AB + bA)„. (8) 

Hi and H2 are respectively the first- and second-order 
terms of the time evolution operator in the coupling con- 
stant g. In case of the large weak value, the approxi- 
mation in Eg. ([6]) is invalid because of the measurement 
backaction [25|. In this Rapid Communication, however, 
we treat the weak measurement as a powerful tool for ad- 
dressing the counterintuitive quantum phenomena seen 
in Refs. [8l-[l2|. Thus, the weak value is predetermined, 
and it corresponds to the so-called fixed value. In this 
case, the measurement backaction can be ignored by set- 
ting a sufficiently small coupling constant. 

First, we calculate the spatial displacement of the 
probe state in the x-y plane. In this case, we employ 
Eq. ^ up to the first order and obtain the spatial dis- 
placements as 

(l)fi - g (Re(i)„ + I Im(S)w) , (9) 
(y)fi = g (Rc{B)^ - I Im(i)„) . (10) 



These results show that the spatial displacements along 
the X and y directions include not only the real part of 
the weak values but also its imaginary part except for 
I = 0. This is because the probe state given in Eq. ^ is 
not factorable in the x and y dimensions so the coupling 
of the observable A (B) with the x (y) dimension of the 
probe also affects its y (x) dimension as the imaginary 
part of the weak value. 

Next, we consider the joint measurement of two po- 
sition observables along the mutually-perpendicular di- 
rections in the probe state. The joint observables in the 
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probe state can be represented by XY and X'^ — Y^. 
The other products of the perpendicular position oper- 
ators can be represented by the hncar superposition of 
these two observables. The expectation values of these 
joint observables can be experimentally obtained from 



the two-dimensional intensity distribution captured by 
the imaging sensor. 

We employ Eq. ([6]) up to the second order in the cou- 
pling constant g and obtain the second-order spatial dis- 
placements as 



,X2_y2 

3' 



P -III -2 



P -III -2 



Rc ({AUB): 



1+2 iAB + BA\ , lA^- 

' ' Re( ) -llml 

4 \ 2 /w \ 2 



B^ 



\{AU^-\{BU?] + 



2 



■Re 



Hm 



AB + BA 



(11) 
(12) 



For the fundamental Gaussian probe state = 0], our 
result is consistent with the result given in Ref. In 
this case, Eqs. ([TT|) and contain only the real part 
of the joint weak values {AB -\- BA)^ and {A^ — B^)w 
For \l\ > 0, however, we obtain additional terms propor- 
tional to their imaginary part as illustrated in Ref. Q 
for the case of |/| = 1 [2^. Thus, the second-order spa- 
tial displacements contain all the information about joint 
weak values. If there are a sufficient number of differ- 
ent outcomes for the first- and second-order spatial dis- 
placements, we can calculate the joint weak values from 
Eqs. dH]) [23|. This feature is not attainable with a 
fimdamcntal Gaussian mode (1 = 0). 



order spatial displacement as follows: 



One useful method for obtaining the joint weak values 
is to take two separate measurements by using two probe 
states with different I values. From a practical point 
of view, the best choice is two probe states with equal 
magnitude but different signs of I since they can be easily 
prepared by using a mirror reflection. In this method, 
the two separate measurements with four types of spatial 
displacements bring eight outcomes, which equals to the 
number of unknown real and imaginary parts of single 
and joint weak values. Thus, it is possible to calculate 
the real and imaginary parts of the joint weak values from 
Eqs. ([9])~ ([T2|) . The advantage of this method is that we 
can use the same interaction Hamiltonian, i.e., the same 
experimental setup; we only need to change the input 
probe state. 



In particular, in case of |^| = 2, the coefficients of the 
first-order weak values (A)^ and (i?)w in Eqs. ((TT|) and 
([T2|) are eliminated. Therefore, only the joint weak values 
remain. By using the two probe states with I = ±2, we 
can extract the joint weak values only from the second- 



Re(ii? + i?i)w = 



{XY)+ + {XY). 



(13) 



lm{AB + BA)^ Y^h - {X^ Y^)- ^ ^^^^ 



Re(A2 - B^) 



Im(A2 - = - 



2g2 



(15) 
(16) 



where ( • )+ and ( • )_ correspond to the spatial displace- 
ments for the I = +2 and —2 cases, respectively. 

Moreover, under the assumption of A^ = B^, the joint 
weak value {A"^ — B'^)vj is also eliminated. Thus, we 
can derive a simple and direct relationship between the 
second-order spatial displacements and the joint weak 
value as 



Re (is -I- BA)^ 
Im(iS + si). 



-(xy)fi, (17) 

9 

'-^{X^^Y^h. (18) 



With these equations, the joint weak value can be di- 
rectly extracted from the spatial displacements by us- 
ing a single LG probe state with \l\ = 2. Although our 
consideration is restricted only to the case of ^ B^ , 
this case includes many experimental setups. For ex- 
ample, the joint weak measurement of Pauli operators 
di {i = x, y, z) on different two-level systems, such as 
A = (Tz (g) / and B = I with the identity operator /, 
is included since they satisfy the property — B^ = I. 

In this Rapid Communication, we derive the full de- 
scription of the weak measurement for the LG probe state 
and provide a simple method for extracting the joint 
weak values only from the spatial displacements of the 
two-dimensional probe state. Our method has a statisti- 
cal advantage over previous ones since information about 
the momentum displacement is not required. Moreover, 
by using the LG probe state with |/| =2, and for the case 
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=13^, we can extract the full joint weak value, i.e. , the 
real and imaginary parts of the joint weak value directly 
from the second-order spatial displacements without any 
additional data processing. 

In our study, we restrict the probe state to the LG 
mode with the radial index p = 0. However, we can use 
other spatial modes, e.g. LG modes with nonzero ra- 
dial index p, Hcrmitc Gauss modes, and hypcrgcomctric- 
Gaussian modes. There remains an interesting problem 
on the relationship between the spatial rotational sym- 
metry of the probe state and the joint weak value since 



the spatial rotational symmetry in our restricted case is 
broken after the weak measurement. Moreover, it would 
be interesting to study the case of probe states given by 
superpositions of different spatial modes. By using such 
a extended probe state, the condition = might be 
relaxed. 
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{AB)w or {]3A)w except for [A, _B] = cl, where c is the 
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